Influence of lattice distortions in classical spin systems 
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We investigate a simple model of a frustrated classical spin chain coupled to adiabatic phonons 
under an external magnetic field. A thorough study of the magnetization properties is carried out 
both numerically and analytically. We show that already a moderate coupling with the lattice can 
stabilize a plateau at 1/3 of the saturation and discuss the deformation of the underlying lattice in 
this phase. We also study the transition to saturation where either a first or second order transition 
can occur, depending on the couplings strength. 

PACS numbers: 75.10 Jm, 75.10 Pq, 75.60 Ej 
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I. INTRODUCTION AND MOTIVATION 

The study of frustrated spin systems continues to be 
a subject of intense research, in particular in low dimen- 
sions where the effect of quantum fluctuations is more 
dramatic, leading to fairly rich phase diagrams. On the 
one hand, one-dimensional frustrated quantum spin sys- 
tems are in general well under control, mainly thanks to 
the availability of powerful techniques like bosonization^ 
and DMRG'^"^. On the other hand, these techniques 
have unfortunately not been successfully generalized to 
the two-dimensional case and there is then a strong need 
for the development of useful techniques to analyze these 
systems^ . 

A standard way to study quantum spin systems is to 
start from the analysis of the classical (large S) limit and 
then try to include the effects of quantum fluctuations in 
a systematic manner^. In certain cases, this procedure 
can lead to a reasonable description of an otherwise in- 
tractable problem. Related to this, the interplay between 
frustration and classical phonons has been shown to lead 
to interesting features, even for the classical spin system 
on the pyrochlore lattice* , like the stabilization of a mag- 
netization plateau at 1/2 of saturation. Then, a natural 
question that arises is whether the classical limit could 
be generally a good starting point to tackle the issue of 
the interplay between frustration and lattice deforma- 
tions and its incidence on the appearance of magnetiza- 
tion plateaux. In the present paper we analyze this point 
by focussing on a one-dimensional Ji — J2 model coupled 
to classical phonons, where both the quantum and clas- 
sical situations can be analyzed and compared. This and 
related problems have been studied long time ago^"^^, 
but to our knowledge the magnetization properties have 
not been analyzed so far. 

The quantum version of this model has been studied 
in a recent article^^, where it has been shown that the 
effects of lattice distortions coupled to a given frustrated 
quantum spin system can lead to new phases, in partic- 
ular to plateaux and jumps in the magnetization curve. 



Although plateaux phases are also present in the pure 
spin system^^, it has been shown that lattice effects can 
lead to the enhancement of these phases under certain 
circumstances. It is worthwhile mentioning that inor- 
ganic compounds like CuGeOa^^ and LiVi205^^'^^ are 
well described by the Ji — J2 model, rendering its study 
both theoretically and experimentally relevant. Values 
for the exchange integrals, such as Ji « IQOK and the 
ratio J2/J1 ~ 0.36, have also been proposed for copper 
germanate^^. 

We shall address the question of whether the effects of 
these lattice deformations can already lead to interesting 
magnetization properties at the classical level. The main 
motivation for the present study is to analyze the origin 
of such plateaux in the particular case of a classical zig- 
zag chain. Although this case is particularly simple and 
the quantum model can be treated using bosonization, 
understanding the role of lattice deformations for clas- 
sical spins could lead to a way to study more involved 
situations, such as two dimensional frustrated systems, 
where analytical techniques arc not so powerful than in 
one dimension as indicated earlier. 

Let us consider the Ji — J2 frustrated chain coupled to 
adiabatic phonons 



-J2 ^ Si ■ Si+2 - ny^ si 



(1) 



In the previous hamiltonian, we chose to modulate only 
the nearest neighbor (NN) interaction term, and to con- 
sider there is no effect on the next to nearest neighbor 
(NNN) coupling. This minimizes the number of param- 
eters in the hamiltonian. We have however checked that 
the inclusion of such a modulation on the NNN couplings 
does not belie our main conclusions. 

In the classical system phonons can be integrated 
out^*, leading to an extra quartic interaction among the 
spins. The effective Hamiltonian, written in units of Ji, 
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reads 

TicS — ^ (^Si ■ Si+i + aSi ■ Si+2 ^ {Si ■ Si+i)^^ 

-hY^St, (2) 

i 

where the fohowing reduced quantities a — J2/ Ji^ Ai — 
-^rj^ and h = H/Ji were defined. Even though one ought 
to study the effect of the elastic constant K and Ai sepa- 
rately, we will focus on the reduced coupling Ai whenever 
possible, reducing the number of parameters to a man- 
ageable size. 

In Section II we study the magnetic phase diagram us- 
ing numerical and analytical techniques. We pinpoint a 
region in the parameter space where a plateau appears 
at Mz — 1/3 only. This should be contrasted with the 
quantum model, which shows in addition a clear Mz = 
plateau in a wide region of the parameter space, and 
another at Mz = 1/2 in a narrower region. Looking 
into the detailed structure of the ground state at these 
plateaux, one can understand this discrepancy in the fol- 
lowing way: the structure at Mz = 1/3 is of the "Up- 
Up-Down" ([/[/D) type, indicating a classical plateau^°, 
while in the Mz = case the singlet structure can be 
identified with a quantum one. 

In Section III we discuss the transition to saturation, 
which is found to be either of first or second order de- 
pending on the ratio between frustration and effective 
lattice coupling. 

II. 1/3 MAGNETIZATION PLATEAU 

Let us analyze the magnetic phase diagram of the 
model (2). In the absence of an external magnetic field 
and when Ai < y/Aa — 1, the ground state is a spiral 
with a pitch angle 9 given by cos6' — l/{Al — 4a). Its 
energy is 

i^spirai = ^ cos - a. (3) 

When Ai > \/Aa — 1 the ground state is Neel ordered. 
The magnetization curves of this system show interesting 
features, which vary depending on the relation between 
a and Ai, as we discuss below. 

In Fig. 1 we represent M{h) for a fixed value of the 
frustration a = 1/2 and different values of the spin- 
phonon coupling Ai . The data was obtained using classi- 
cal Monte-Carlo (MC) based on the usual Metropolis al- 
gorithm. Starting at high-temperatures we perform sev- 
eral thousands of MC sweeps, and then cool down the 
system to a fraction of the initial temperature. This pro- 
cedure is then repeated, slowly annealing the system to 
zero temperature. We observe that a steady magnetiza- 
tion plateau at 1/3 appears as soon as the coupling to 
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FIG. 1: Magnetization curves M{h) for A'^ = 30 spins with 
a = 0.5 and Ai = 0.0, • • ■ , 0.8 in steps of 0.2. Periodic bound- 
ary conditions are applied on the chain. The system is grad- 
ually cooled to T ~ over 3 x lO" Monte-Carlo sweeps (color 
version online). 

the lattice is slightly turned on, whose length increases 
with Ai. One can notice that the way the system enters 
the plateau from the low-field side and eventually satu- 
rates differs depending on the effective lattice coupling 
Ai. For Ai > 0.6 the two are first order transitions. An- 
other interesting characteristic seen in Fig. 1, is that all 
curves represented (except one) cross at the same field 
/ix « 3.35 for which Mx = M(/ix) « 0.745. We shaU 
discuss this point at the end of section III. This brief 
overview suggests that the coupling with the phonons 
stabilizes the state at Mz = 1/3. Since the plateaux are 
observed at zero temperature, we can fairly assume that 
this effect is energy driven. 

After this numerical preamble, it is time to derive some 
analytical predictions on the characteristics of the mag- 
netization plateau. For this purpose, we need to find out 
which states describe the system in the low and high-field 
regions around Mz = 1/3. It is not surprising that in the 
plateau phase the system adopts the UUD state, but it 
will be a crucial point in our discussion as we shall see 
later. In this state, the spins are aligned along the z- 
axis, two up spins alternating with one down spin which 
is precisely the structure seen at the 1/3 plateau in the 
quantum model^^. 

The classical MC data indicates that the situation in 
the low-field region is more complicated. On the one 
hand, the transition to the UUD state can occur at a 
very low field, where the system is not far from its zero 
field ground state. Then, there is no small unit cell struc- 
ture providing a good description of the system, as the 
spiral structure still prevails. On the other hand, when 
the transition is smooth in the low-field region, a plau- 
sible assumption is to consider that the system adopts 
a coplanar "Y" configuration parameterized by a single 
angular degree of freedom 6 (see Fig. 2). The unit cell 
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energy for this state reads 

Ey{9) = (l + a)(2cos6'(cos6i-l)-l) 

~h{2cose~ 1). 



(4) 



This expression can be minimized for any set of the pa- 
rameters h,a and Ai. As the magnetic field increases, 
the solution will eventually yield 9 = corresponding to 
the UUD state. This configuration is always a solution 
of deEyiO) = 0, but it is only a minimum of the energy 
when h>hY = l + a — 3Ai^ . We should emphasize that 
this discussion only makes sense whenever hy is positive. 
For a given value of the magnetic field, there can be other 
solutions satisfying 
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X2. 

(5) 

where X = sin 9 (assuming 0036* > 0). The study of Eq. 
(5) boils down to finding the sign of a polynomial expres- 
sion. Introducing A = 2(1 -fa) — llAi^, we can show 
that when A > there is exactly one more extremum 
of the energy for h < hy and that it is always a mini- 
mum. This solution becomes precisely the UUD state at 
h = hy. Under these assumptions, we can conclude that 
the critical field for which we recover = 1/3 from the 
low field regime is 



hci = l + a-3Ai^, A>0. 



(6) 



This can be compared to our MC results. For instance, 
the data for a = 0.5 and Ai = 0.4 (solid pink curve in 
Fig. 1 ) allows us to obtain a precise estimate for hd at 
T w 0. We get h^i = 1.02 ± 0.01. For this set of pa- 
rameters, A is positive so that we are ruled by the pre- 
vious assumptions. The analytical expression (6) yields 
hci = 1.02, which is in excellent agreement with the sim- 
ulations. 

For A < there can be up to two extra solutions 
when h > hy. As there is always one solution that 
never turns out to become UUD for a certain value of 
the magnetic field, we ought to perform a detailed com- 
parison of the two solutions' energies in order to con- 
clude. We shall not step further into this discussion, 
which can nevertheless be conducted numerically using 
the previous analytical expressions. For instance we per- 
formed it when a = 1/2, Ai = 0.6, leading to hd ~ 0.46. 
This is in good agreement with the MC data which gives 
hci ~ 0.47 ± 0.01 (dotted curve in Fig.l ). It can be un- 
derstood from the previous discussion that hy is always 
a lower boundary of the critical field: 



hci > 1 



3Ai^, A<0. 



(7) 



If we increase Ai while keeping a fixed, hy eventually 
becomes negative (as it is the case for Ai = 0.8, dash- 
dotted blue curve in Fig. 1) and we can generally not 
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FIG. 2: Configurations observed in the low-field, the 1/3 mag- 
netization plateau and high-field regions. The chain is viewed 
in the ly-plane. The arrows denote the projection of Si in 
this plane, whereas the circles represent the Sz component 
(red for Sz > 0, yellow otherwise and radius proportional to 
\Sz\, color version online). The parametrization of the states 
is given for each configurations. 




FIG. 3: Magnetization curves M{h) for A'' = 30 spins with 
Ai = 0.3 and a = 0.0, 0.4, 0.6, 0.8. Periodic boundary condi- 
tions are applied on the chain. The system is gradually cooled 
toT over 3 x lO** Monte-Carlo sweeps (color version on- 
line) . 



conclude using this small unit cell configuration. The 
reader should keep in mind that the regime where Ai 
becomes large is not well described by our initial hamil- 
tonian (1) since in that case one should include the effects 
of the lattice also in the NNN interactions. 

We shall now focus on the state observed in the high- 
field region to find the corresponding upper critical field 
hc2 above which the plateau disappears. By comparing 
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hci to hc2, we should be able to conclude on the existence 
of the 1/3 magnetization plateau. In the upper critical 
region, the situation is far more under control. The sys- 
tem can be seen to be well described by a 3-spins coplanar 
"canted" configuration with two degrees of freedom (see 
Fig. 2). The energy of such a configuration is given by 

-^(l + 2[/(^?i,e,)2) 
-/i(2cos6'i - COS6I2), (8) 

where U {91,62) = sin 6'i sin 6*2 + cos 9i cos 02 ■ 

The configuration UUD, which corresponds to 9i = 
O2 = 0, is always a critical point of the function 
£'cantod(^i, 6*2)- A closcr look at the second order deriva- 
tives with respect to di and 62 shows that it is a local 
minimum only for < h < 1 + a + Ai'^. The other critical 
points satisfy the following set of equations 

Y ^ 2X, (9) 

h = (l + a + Al (2X^ + Vl - ^Vl - 4X2(71(72)) 

X (2y/l-X^ - y/l - 4X2(71(72) , (10) 

where 

sin6li = X,cos6'i = aiy/{l~ X^), (11) 

sin 02 = r,cos02 = ^2^(1-^')- (12) 

The quantities (7i,(72 = ±1 account for all the possible 
signs of both cosines. We see from Eq. (9) that there is 
a strong constraint on {61,02) verified regardless of the 
values of the couplings. At h = 1 + o; + Ai^ , Eq. (10) 
admits only one solution which turns out to be UUD. 
For larger value oi h, UUD can no longer be a critical 
point, which implies 

/ic2 = l + a + Ai^ (13) 

corresponding to the exit of the plateau in the high-field 
region. From this discussion it can be concluded that 
whenever our assumptions are correct, there is a plateau 
at Mz — 1/3 of length A/11/3 — AA\ starting at hd- 
This result has been checked to be consistent with the 
MC computations and the analytical value of hc2 matches 
the value estimated from all the curves in Fig. 1. 

There is one more question we need to address: for 
which set of parameters {Ai , a) can we observe this 
plateau? Under the previous assumptions regarding the 
states observed in the low and high-field regions, we can 
conclude it exists for any Ai > 0. Yet the system can 
not be described in such a manner for all values of a and 
Ai. Working at a fixed lattice coupling Ai — 0.3, we 
were able to obtain some magnetization curves varying 
the frustration a. Some of those curves are plotted in 
Fig. 3, which clearly shows that there is only a narrow 
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FIG. 4: Qualitative {A\,q) piiase diagram (color version 
online). The filled area between the two full curves cor- 
responds to the region of the parameters space where the 
Mz = 1/3 plateau is observed. The hatched region between 
the two dashed curves is the region of the parameter space for 
which our approach is no longer observed to be fully valid. 
Also represented the limit between Neel and Spiral ground 
states at /i = (blue dash-dotted line) and the region where 
1 -|- a — < (dotted orange line). 



region in a where the plateau is observed. A precise an- 
swer to the previous question is rather challenging, and 
we shall first try to discuss this point in a more qualita- 
tive manner before adopting a more precise strategy. At 
= 1 /3, we can of course expect to see a lot of different 
configurations, depending on the values of the couplings. 
However, the MC simulations suggest that the plateau al- 
ways corresponds to the UUD configuration. This state 
is perfectly collinear, minimizing the quartic contribution 
to the effective hamiltonian (2). For instance, it can be 
seen numerically that for ta; = 1/2 with no coupling to 
the lattice, the system reaches 1/3 magnetization in the 
UUD configuration. Even a small positive value of A\ 
will then stabilize the UUD state enough for it to be sta- 
ble when the field is slightly increased. On the opposite, if 
one antiferromagnetic coupling dominates the other, the 
system will be in a different state at = 1/3. In the 
extreme case where o: s:! for instance, the system will fa- 
vor Neel order in the xy plane, each spin having the same 
2;-axis projection = 1/3. This layout already trades 
off some collinearity in favor of magnetic field alignment. 
There is no surprise that this trade-off will be further 
enhanced as the magnetic field is increased, so that no 
plateau should be observed. 

A more accurate way to tackle this issue is to start 
from the /i = spiral ground state and ponder over the 
state adopted by the system when the magnetic field in- 
creases. We have already performed part of this task ear- 
lier, suggesting that the system slowly moves to a "Y" 
configuration, whose out of plane components make it a 
"precursor" of the UUD configuration. Another plausi- 
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FIG. 5; Magnetization curves M{h) for A'' = 30 spins with 
Ai = 0.3 and a — 0.4 around = 1/3 for different temper- 
atures T = 0, 0.04, 0.08 in units of Ji. Averages are computed 
on 2^^ sweeps through the lattice after an initial 2^* sweeps of 
thermalization. Increasing the temperature quickly destroys 
the plateau observed at Mz = 1/3. 



ble solution is that the spins, while keeping their spiral 
structure in the xy-plane, all acquire the same Sz projec- 
tion. In this case, the n}^ spin reads 

Sn = (Vl - cosine), Vl-z^ sm{ne),z^ , (14) 

with cos 6* = l/{Al — ia). The energy per site. 



(15) 



4(z^+(l-z^), 



hz. 



can be minimized with respect to z to find the lowest 
energy configuration at a given magnetic field. Our idea 
is to perform this minimization ai h — hd , and to see if 
the corresponding configuration is of lower energy than 
UUD at the same field. If so, the system will not enter the 
plateau at hd , and of course as mentioned the previous 
paragraph no plateau should be observed. For a fixed 
value of Ai , we can determine the range in a leading to 
UUD at hci- The roots of the polynomial equation are 
evaluated numerically, from which we sketch the phase 
diagram represented in Fig. 4. This approach only makes 
sense when we have a precise value for hd , which we saw 
is the case if Ai is not too large {Ai < 0.5 from the MC 
data). We notice that the diagram is in agreement with 
the situation depicted in Fig. 3, as well as the one in 
Fig. 1 when Ai is not too large. The most remarkable 
feature is that for an arbitrary small yet strictly positive 
Ai , one can find a value of a for which the plateau phase 
is observed. 

The effect of temperature on the magnetization plateau 
is potentially important as an "order by disorder" effect^^ 
could further stabilize the plateau. We investigated this 



point by performing our MC simulations at different tem- 
peratures, without annealing the system. A sample is 
given in Fig. 5, and in general we observed no remarkable 
features. The increasing thermal fluctuations quickly de- 
stroy the plateau. We should also mention that we ob- 
served no strong finite size effects in the numerical sim- 
ulations, which is why we were always able to work on 
systems with less than a hundred spins. 

We conclude this section by focusing on the lattice de- 
formations. Until now, we studied the effective spin-only 
hamiltonian (2) which embodies the straightforward ana- 
lytical approach to the problem. It is however important 
to get more insight on the structure of the lattice defor- 
mation inside the plateau phase. For that matter, we 
modified our MC algorithm to take into account the lat- 
tice degrees of freedom as well. Starting from the hamil- 
tonian (1), we used the Metropolis algorithm for both the 
spin positions and orientations, applying periodic bound- 
ary conditions on the chain. We studied the normalized 
histograms of the displacements Si at finite temperature. 
We fixed a = 0.4 and Ai = 0.3, the same values used 
in Fig. 5 to allow a direct comparison between the two 
figures, and selected the magnetic field so that the sys- 
tem is at Mz ~ 1/3. Besides the value of Ai = Ai/\/k, 
we need to give K, the spring constant in (1), a sensi- 
ble value. We took K = 10'^ Ji, large enough to make 
sure the displacements remain small. This corresponds 
to Ai Ri 9.5. We mention that both K and Ai are of the 
same order of magnitude as the one for a more complex 
two dimensional material such as SrCu2(B03)2^ and that 
they can be considered at least as "realistic" for copper 
germanate or lithium vanadate^^. The results are given 
in Fig. 6. 

We see that the lattice deformations are not uni- 
form and that their histogram presents two peaks at 
T = O.OlJi. They are centered around a negative and 
positive value of the displacement Si. This suggests that 
the underlying deformation consists of "UDU" trimers 
(Up-Down- Up) on the chain. Let us introduce S-^- the dis- 
placement between two consecutive trimers and S- the 
displacement between the down spin and its two nearest 
neighbors inside the trimer. The energy of this unit cell 
is given by: 



E = SKSl + iJiAiS^ -J1 + J2-H, 



(16) 



where the periodic boundary conditions imply = 
—2S-. Minimizing the energy, the deformation should 
become 



-2S^ 



AAiJi 
3K 



(17) 



at T = OK. Going back to Fig. 6, at T = O.OlJi the 
distribution clearly exhibits two peaks and we can see 
that they are almost centered around S^ and S- respec- 
tively. The ratio between the height of the two peaks 
is about 2, a consequence of the fact there are twice 



6 



— r = o.oiJi 




-0.04 -0.02 0.02 0.04 

FIG. 6: Normalized lattice displacement histograms for A'' = 
30 spins with a = 0.5, Ai = 0.4 {K = 10^ Ji) and h = 1.5 at 
different temperatures. 300 points were used in the interval 
[—0.1,0.1]. The two dotted vertical lines correspond to the 
T = limit calculated in the text. The data was obtained 
using a direct classical Monte-Carlo for the hamiltonian (1). 

as many up spins than down spins in the UUD state. 
Those results seem to validate the trimer scenario at 
low temperature. When the temperature increases to 
T = 0.04Ji, the peaks start to overlap, betraying the 
gradual destruction of the plateau already seen in Fig. 6. 
Finally at T = O.O8J1, we end up with a single peaked, 
almost gaussian, distribution: the plateau eventually dis- 
appeared. We end up by stating that expectation value of 
the displacement is always zero as the periodic boundary 
conditions applied ensure the length of the chain remains 
fixed throughout the simulation. 

III. TRANSITION TO SATURATION 

The study of the upper critical magnetic field yields 
another interesting result: we can get a precise picture of 
how the system eventually reaches saturation. This re- 
sult can be foreseen using classical MC, which shows that 
the canted state describes the system quite well even for 
h > hc2- A close look at Fig. 1 shows that two dif- 
ferent behaviors of the magnetization between hc2 and 
the saturation value are observed. For different values of 
spin-phonon coupling, the system can undergo a first or 
second order transition to reach saturation. We are go- 
ing to demonstrate that this result can be derived from 
energetic considerations on the canted state. From Eq. 
(8), we see that the saturated state, reached for 61 = 
and 02 = TT, minimizes the energy for a magnetic field 
greater than hc3 — 3{1 + a — Ai^). This imposes a lower 
boundary on the saturation field hu- We assume that 
the couplings Ai and a are such that hc2 < hcs, a sit- 
uation where the previous discussion on the existence of 
the 1/3 magnetization plateau still holds. To be consis- 



tent with the state of system for h > hc2, we set cti = 1 
and let a2 = —u take the values ±1 so as to be able to 
move from UUD to saturation continuously. At a given 
magnetic field, one can obtain the corresponding critical 
configurations by finding the roots of Eq. (10). This task 
reduces to the study of the two functions ha- : 

ha{X) = (l + a + Aj (2X^ - - ^ Vl - ) 
X (2\/l-X^ + \/l - AX^aj . (18) 

Their roots can be determined graphically for a fixed 
field /ig as they are the values of X for which the line 
h = /iQ intersects h^iX). The "low magnetization" func- 
tion will give us solutions with one spin still point- 
ing down, whereas the "high magnetization" function 
will give us states where all the spins have a positive Sz 
component. Fig. 8 is a plot of both functions for two 
sets of values a, Ai. In both cases, the curves for 
and /i„ join at hc4 = \/3(l + a + Ai/2) (colored dots 
in Fig. 8). For this value of the magnetic field, the root 
of ha- corresponds to a configuration in which one of the 
three spins lies precisely in the xy-plane. Two possible 
behaviors are observed. For instance when a — Ai = 0.5 
(green curves in Fig. 8), we see that for a fixed mag- 
netic field h e [hc2, hcS\ there is only one critical point of 
the energy, which can be shown to be a minimum. We 
are able to follow easily the state of system as the mag- 
netic field increases. The two up spins first slightly tilt 
to let the down spin reach the xy-plane and then they all 
progressively align along the z-axis while still satisfying 
Eq. (9). The three-spins unit cell configuration smoothly 
goes from UUD to saturation. 

For Ai — 0.8 (black curves in Fig. 8), the "high mag- 
netization" function /i^ (plain line) presents a maximum. 
In this case, three states are potentially competing for h 
between hcs and its maximum value: the saturated state 
and the two roots of We ought to compare their 
energies to conclude, but it is not surprising that the 
outcome can be a first order transition to saturation. We 
numerically solved the analytical equations involved to 
get the magnetization curve from the exit of the plateau 
to saturation for Ai — 0.8. The saturation field we ob- 
tain is hu « 3.17063, for which the system jumps from 
Mz ~ 0.58434 to saturation. The comparison between 
this minimization and the MC data is given in Fig. 7 
and shows the excellent agreement achieved. 

A more in-depth study of the ha functions' extrema 
allows to work out the range in {a,Ai) for which the 
transition to saturation is of first or of second order. The 
former, which are related to the existence of a non trivial 
maximum in h-^^, occur only if 1 < {1 + a) /A^ < ^■ in 
agreement with our numerical observations. 

Finally, it can be pointed out from Eq. (18) that for 
X — 1/ \/5 the function h^ does no longer depend on the 
coupling Ai. This state (if reached) will be the minimum 



7 



0.8 




X )( X X 



X M.C. 
— Theo. 



2.4 



2.6 



2.8 



3.2 



FIG. 7: Upper region of the magnetization curve for a — 0.5 
and Ai = 0.8 at T = 0. The full line represents the results 
obtained using the analytical minimization of the "canted" 
state energy, crosses the data obtained by Monte-Carlo. 




FIG. 8: "Low magnetization" (dashed line) and "high mag- 
netization" (plain line) curves as a functions of X = sin 9i for 
a = 0.5, Ai = 0.5 (green) and a = 0.5, Ai = 0.8 (black). 
The typical shape of the minima's unit cell along the green 
curves are depicted in the different magnetic field regions. 
The intersections (crosses) with the line h — ho (red) gives 
the competing critical point at this magnetic field. The col- 
ored dots correspond to the minimal configuration with one 
spin in the xy-plane. 



of the canted configuration energy for a magnetic field 



(19) 



At this field the magnetization is = = 
5/{3^/5). This explains why for our selection of param- 
eters, all the curves except one in Fig. 1 cross at a field 
whose estimate, given in Section II, coincides with hx- 
Regardless the value of Ai , if the system is not saturated 
at hx then its magnetization will always be My. 



IV. SUMMARY AND CONCLUSION 

The effect of lattice deformations at the classical level 
in a frustrated spin system has been illustrated work- 
ing on a simple Ji — J2 spin chain coupled to adiabatic 
phonons. We provide an overall picture of the magneti- 
zation properties for a large set of the parameters a, Ai 
introduced in our model. We have found that a plateau 
at AIz — 1/3 is present in certain region of the parameters 
space, while no other plateaux are observed. Frustration 
is a necessary ingredient, as the plateaux can only arise 
when the zero field ground state is a spiral. The other 
ingredient, the coupling to lattice deformations, is such 
that for an arbitrary small yet strictly positive Ai, one 
can find a value of a for which the plateau phase is be- 
held. Further increasing Ai will broaden the region in 
the parameter space for which the plateau occurs, until 
the effective coupling is no longer mild enough for our 
analytical approach to be valid, even if a numerical ap- 
proach is still achievable. It should be emphasized that 
the stabilization mechanism is purely energy driven and 
triggered by the quartic interaction induced by the lattice 
coupling in the effective hamiltonian (2). The underlying 
lattice deformation shows the chain is made of "UDU" 
trimers inside the plateau phase. 

The absence of plateaux at — and = 1/2 in 
the classical model as compared to the quantum case can 
be understood by analyzing the ground state structure of 
the plateaux in the quantum case. It is only for Mz = 1/3 
that one observes a classical type of spin configuration, 
of the UU D type, while in the other cases a quantum 
state is apparent. 
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